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Методика ошнки оптимальних обсягв 
статистичних випробувань байесовським методом 


Одним из основных требований к испытаниям исследуемой аппаратуры, повышающим их эффективность, 
является уменьшение числа испытаний при заданной достоверности статистической информации. На основе 
байесовского метода, позволяющего учесть априорные данные статистических исследований, составлен 
алгоритм оценок результатов статистических испытаний, обеспечивающий заданную достоверность с 


менышим числом испытаний. 
Ключевые слова: достоверность информации, оптимизация, эффективность. 
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Введение 


При испытаниях аппаратуры возникает задача достоверной оценки её параметров, 
достигаемой с минимальным числом испытаний. Байесовские методы, применяемые в 
статистике, основаны на учете априорной информации, вследствие этого позволяют 
существенно сократить число статистических испытаний при заданных требованиях на 
достоверность статистических оценок. 

На основе сопряженных распределений вероятностей можно составить алгоритмы, 
которые снижают трудоемкость вычислений и позволяют сравнительно простым спо- 
собом оценить основные характеристики аппаратуры по данным статистических испы- 
тании. 
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Для определения оптимального числа испытаний в качестве основного критерия 
оценок исследуемых параметров применяется байесовский риск, который при квадратич- 
ной функции потерь, учитывающей возможность отклонения случайных значений ис- 
следуемого параметра от его среднего значения, определяет суммарную дисперсию ис- 
следуемой, случайной величины, включаемую как дисперсию самого математического 
ожидания относительно его истинного среднего значения, так и дисперсию случайной 
величины относительно математического ожидания, а также изменения дисперсии, за- 
висимые от числа испытаний. При этом важно учесть также корреляционную связь слу- 
чайных величин математического ожидания и меры точности исследуемого распреде- 
ления вероятностей, определяемой величиной, обратной дисперсии распределения. 

В процессе минимизации байесовского риска можно получить оптимальные объемы 
требуемых испытаний на двух этапах проведения таких испытаний, и взаимосвязь этих 
величин можно оценить с учетом влияния ряда факторов, оказывающих воздействие на 
объемы испытаний на первом и втором этапах их проведения. 


Основная часть 


Согласно теореме Байеса [1], [2] условная функция плотности вероятностей 
(ФПВ) Р( / х) для вектора параметров 0, который включает математическое ожи- 


дание и дисперсию, считающиеся также случайными величинами при условии полу- 
чения вектора Х случайных наблюдений, может быть выражена следующей зависи- 
мостью: 


Р@).Р(Х /0) 
[Р®). Р(х/о 40° 


Ре /х)= 


где Р®) и Р(Х / 0) — соответственно ФПВ распределения вектора параметров 


и условная ФПВ распределения вектора х наблюдаемых значений при заданном зна- 
чении вектора 9 параметров, которая называется функцией правдоподобия. 
Считается, что Р(х/ 0) есть априорное условное распределение вероятностей, а 


Р( / х) — апостериорное условное распределение. 


Планирование и анализ эксперимента на основании теоремы Байеса существенно 
упрощается при введении сопряженного стандартного семейства распределений вероят- 
ностей вектора 9 параметров, обладающего следующим свойством [2], [3]. 

Если априорное распределение 9 принадлежит этому семейству, то при любом 
объеме выборки п и любых значениях наблюдений в выборке апостериорное рас- 
пределение вектора 9 параметров должно также принадлежать этому семейству. 

Таким образом, сопряженное семейство распределений предполагает особую 
связь, которая должна существовать между семейством распределений вероятностей 
параметров и семейством распределений вероятностей наблюдений. 

На основании изложенного в работах [2], [3] предлагается ряд способов оценки 
параметров в зависимости от вида распределений вероятностей при получении пов- 
торных выборок наблюдаемых значений. 

Рассмотрим одномерное нормальное распределение случайной величины при 
неизвестных параметрах: математическом ожидании и дисперсии. С целью упрощения 
математических преобразований удобнее в качестве меры рассеивания случайной вели- 


чины вместо дисперсии с - нормального распределения рассматривать меру точности 
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т нормального закона, определяемую как величину, обратную дисперсии т =1/с 
Полагая, что хи..... х„ — повторная выборка из нормального распределения с неизвест- 


ными значениями среднего М и меры точности К ‚ а априорное совместное распреде- 
ление М и К определяются в виде: 
Условное распределение М при К = "(к > 0) — нормальное со средним ци и мерой 


точности т ‚ причем — < < ци <с0л > 0; маргинальное распределение К есть гамма- 
распределение с параметрами © и В, где а>0 и В >0, которое определяется 
выражением для вероятностной плотности: 


69) 9 
(")= Ге)" з е реа, › 
0 в остальных случаях 
где математическое ожидание Е(ь) случайной величины Г и соответственно 
дисперсия Б(+) равны: 


и а. 


согласно [2], [3] апостериорное совместное распределение М и К при х=х, 
г = 1,2,....п определяется следующим образом: 
— условное распределение М при К =г есть нормальное со средним и’и 
мерой точности т ' = ( + пу, где 
‚ учи 


НЕТ, 
т-п 


— маргинальное распределение К есть гамма-распределение с параметрами 


} п 
о =а+— 
2 


в'=В+5 (= ри (1 
2(° + п) 
Полученное сопряженное семейство распределений вероятностей обладает ря- 
дом важных особенностей: 
| - априорное маргинальное распределение М есть распределение Стьюдента 
ох 


2 — апостериорное маргинальное распределение Стьюдента получается заменой 
их, 9 и В их апостериорными значениями, причем: 

— число степеней свободы 29 +п апостериорного маргинального распреде- 
ления М не зависит от наблюдаемых значении; 
— параметр сдвига п и точность т зависят от наблюдаемых значений, 

3 — параметры М ик всегда зависимы; 

4 — если значение меры точности К распределения вероятностей велико, то 
наблюдение дает очень точную информацию о значении среднего М, если же 
значение К мало, то одно наблюдение дает мало информации о значении М. 

При условии, что исследуемые параметры системы имеют нормальное распре- 
деление, можно на основании приведенных выше соотношений составить алгоритм по- 
следовательных оценок средних значений ЕМ | ЕК), а также дисперсий Ым), В), 
математического ожидания М и меры точности К. 
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В качестве априорных исходных данных Ё (м ), Е (к) р БМ ) И БК) можно ВЗЯТЬ 


приближенные значения, которые с увеличением числа последовательных наблюде- 
ний не будут оказывать влияние на конечные результаты оценок. Вычислительные 
процедуры алгоритма можно представить по шагам. 


Шаг 1. По исходным данным Е (к) и Б(К) необходимо определить параметры 
гамма-распределения согласно выражениям: 


Ее в 8) 


Б(в) ^^ обв). — 


Шаг 2. Поскольку случайная величина М имеет распределение Стьюдента со 


ю о .т в 
20 степенями свободы и мерои точности р. ‚ ТО из выражении (2) и соотношения 


[2]: р(м )= Вх п. (@ 1) следует 


ы Е(К) 
[е(®)- р®рм) (3) 
и=Е(М). м 


Шаг 3. По данному нормальному распределению производится повторная вы- 
борка из п наблюдений и определяются: 


х=ТУх, => (хх). 
И = 1=1 


Шаг 4. Определяются параметры ц’л'’,о’ и В’ апостериорных распределений 


М и К согласно следующим выражениям: 
—1 


и’ = (и +их)( +п) д’ЕТ+П 
Иан”, "В+ (5) 


Шаг 5. Определяются апостериорные значения математических ожиданий и 
дисперсий оцениваемых параметров М и К: 


ЕАМу-и" ВМ Ц, (©) 
Е(в=® Б\ив)=“ (1) 


В ‚ (В ‚? 

Теперь можно произвести анализ полученных результатов, например, способом 
интервальных оценок. 

Интервальная оценка среднего значения Е (м ) может быть определена по табли- 
цам распределения Стьюдента при следующих исходных данных: 

— число степеней свободы равно К = 24 ; 

— параметр сдвига равен ц; 


— мера точности определяется выражением ый: . 


В качестве значений ©,В, и ит для интервальной оценки М берутся последние 
расчетные данные, полученные согласно описываемому алгоритму. 

Задаваясь, например, доверительной вероятностью у, можно по таблицам распре- 
деления Стьюдента и полученным исходным данным определить доверительный интер- 
вал ( №, № ‚) ‚ в котором находится полученное среднее значение М согласно выражению: 
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ох |? 
Р/-в5| | (М-н)<а|=у: (8) 


В 


Условие (8) можно записать следующим образом: 


Р,‚[М№, М < М, =ъ, (9) 
1 С 
2 2 
ре М. =и-Ь ее (10) 
ОТ ОТ 


причем :( Е ) — квантиль, # — распределения Стьюдента. 
ох 
Выражение (9) означает, что для данного интервала, выбранного из совокупности 
интервалов, вероятность содержать среднее значение М равна у. 
Из анализа равенств (10) следует, что с увеличением числа выборок величина 
доверительного интервала уменьшается. Если задаться требуемым значением довери- 
тельного интервала ( №, № И полученного при заданной доверительной вероятности у, 


то можно повторно производить байесовские оценки параметров согласно изложенному 
алгоритму по опытным данным до тех пор, пока, например, не будет получено заданное 
значение доверительного интервала при заданной доверительной вероятности. 

Разумеется, можно применять и другие способы оценивания значений пара- 
метров, полученных согласно данному алгоритму. Важно только, чтобы было вве- 
дено определенное решающее правило, по которому можно было бы установить: 
продолжать ли последовательные статистические оценки параметров по данному 
алгоритму или же их больше не производить. 

Таким образом, если полученные статистические значения параметров и их ин- 
тервальные оценки являются сомнительными, то процедуры статистических оценок по 
данному алгоритму следует повторять. Для этого производится еще и наблюдений 
из данного нормального распределения и повторяются шаги 2 — 5 описываемого алгорит- 
ма, но вместо исходных априорных значений параметров ЕМ ) Бм ) ЕК), ГХК), ЦТ ‚0 
и В необходимо применять апостериорные значения Е(М ) р(м ) Е(В) (Е), ит’, о 
и В’, полученные согласно выражениям (5) и (6). 

Анализ математических зависимостей позволяет сделать следующие выводы: 
исходные априорные данные по Ё (м ) р(м ) Е (к), (к) не влияют на окончательные 
значения оценок; данные измерений, последовательно используемые при оценках, 
соответственно снижают дисперсии Б(м ) И (К) математического ожидания Ми 
меры точности К ; среднее значение Е (в) меры точности зависит от выборочной дис- 
персии 5 каждой серии измерений; доверительные интервалы оценок среднего значе- 
ния, Ё (м ), полученные при доверительной вероятности у = 0,95 , существенно умень- 
шаются с увеличением числа серий измерений 1. 

Таким образом, предложенный алгоритм оценки параметров позволяет производить 
последовательные байесовские статистические оценки таких параметров, как математи- 
ческое ожидание и дисперсия при нормальном законе распределения вероятностей. Такой 
способ позволяет существенно сократить необходимое число испытаний и контроля ап- 
паратуры. Предложенным алгоритмом можно пользоваться и при двухмерном или трех- 
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мерном нормальном законе распределений вероятностей, когда наблюдаемые статисти- 
ческие данные не коррелированны. Если исследуемые статистические данные коррели- 
рованны, то можно предварительно произвести их декоррелирование способом поворота 
системы координат [4]. 

Как уже указывалось, задача заключается в том, чтобы с помощью статистических 
испытаний байесовским методом оценить векторный параметр, включающий в качестве 
компонент неизвестные значения математического ожидания и дисперсии исследуемого 
нормального распределения. При этом совместное распределение случайных величин ма- 
тематического ожидания и дисперсии есть нормальное гамма-распределение. 

При оценке векторного параметра принимается матричная квадратичная функция 
потерь вида [2]: 


П@,а)= @ а) дб -а) 
где А- симметрическая неотрицательно определенная матрица. 


В качестве критерия оптимальности принимается минимум байесовского сред- 


него риска. Оптимальная по этому критерию оценка вектора 0 есть вектор средних 
значений маргинальных апостериорных распределений математического ожидания и 
меры точности. Байесовский риск при заданном априорном распределении после 
проведенных и испытаний равен матрице, определяемой произведением матрицы А 
на ковариационную матрицу [2] (т, г/ х) : 


Р(иьт.)= {АЖ (т,г/х}}, 
Р(т/х) К(т,г/х 
К(т,г/х) Ох) 


ческого ожидания т и меры точности г соответственно; К(т,г"/х) — ковариация 


где \(т,"/х)= ‚ Б(т/х), Б("/х) — дисперсии математи- 


этих же параметров для апостериорного условного распределения. 
Элементы матрицы (т, г/ х) представляют средние значения дисперсий и КО- 


вариаций параметров ти г. Полагая, что элементы матрицы А равны Ти учитывая 
при этом, что потери, связанные с увеличением объемов испытаний ий; ; возрастают 


пропорционально коэффициентам с,, можно получить следующее выражение для 
общего байесовского риска: 


Р(в, ть} = В{т/х)+ГХ"/ Х)+2К(т,г/Х)-+ ст + сть. (11) 
С увеличением объемов испытаний И, и И, дисперсии, входящие в формулу 


(11), имеющие априорные большие значения, будут уменьшаться, а потери, обуслов- 
ленные (<, будут возрастать. Таким образом, скалярная функция р (и, п.) имеет ми- 


нимум, и оптимальные объемы испытаний п ,П› определяются минимально возможным 


значением риска (11). Расчеты показывают, что последним слагаемым в выражении (1) 
можно пренебречь. Тогда на основании алгоритма, после несложных преобразований 


получим простые рекуррентные зависимости для параметров 7 ,%,, р и И, после К 


серий испытаний с объемами выборок И, : 


К 1 « 14 
тети 059+ Хи В. В+ Утя, и» 
пи 11 


11 
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Га —1 
(иж . (13) 
ы 1=1 
На основании (6) и (12) легко получить рекуррентные выражения для дисперсий 
р, (т) и р. (+) после 2-х серий испытаний с объемами выборок пи и»: 

р, (т/х )= В + п: 5: + 53) (24 -2)+ МС +2а -2)+ М 2", (14) 
р,("/х)= 2@а +М)(2В +п151+ 1252) › (15) 
где М=и, +п,, а параметры а,В и т определяются выражениями (2). Таким 
образом, выражения (14) и (15) позволят достаточно точно оценить условные апосте- 
риорные дисперсии математического ожидания т и меры точности г исследуемого 

нормального распределения случай величин Х после проведенных испытаний. 
Как уже указывалось, условная апостериорная плотность вероятностей не- 


известных значений математического ожидания т и меры точности г нормального 
распределения определяется законом нормального гамма распределения: 


(и, = Ви Г@ "т Иа - в оный И (16) 


где Г@) — гамма-функция. 
Для заданного закона распределения вероятностей (16) корреляционный момент 
К(т,"/х) определяется известной зависимостью : 


К(т, г/х)= = | т- «в онл (17) 
—00 
Интеграл (16) легко раскладывается на 4 интеграла: 
К(т,г!х)= п- Г. 15+ 1, (18) 
где |= [Г рио(т, Г оатаь им 
—® 0 —о 0 


= | [ть (т, " / х)атаг а в | [У(т, г/ х)атаг. 
о 0 20 
На основании (16) интеграл [, можно записать следующим образом: 


Г —УЙ [ те) е-^ ата! , 


—® 0 
ее с 1 ы а 
где а В [= Вт (2 2 Г(а}' (19) 
Интеграл ], является табличным [5] и его можно представить таким образом: 
в Э е 
Е] Це а (20) 


подставляя (19) в (20) и производя необходимые преобразования, нетрудно получить: 


та мо 
В] п?-этичт? ав | тат › (21) 
о т 


где =2Е. Г + 2] ) Е (22) 


т 
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Интеграл (21) — табличный [5] и может быть записан в виде: 


—1 

[= [и В(2;2а +1). „Е Ба + ‚за ое в [+2422 (23) 
т 
где В(2;2а + 1) — бета-функция, определяемая выражением: 
В(2;2а +1) = Г(2)- Г(2а, +1) Га +3) (24) 
—1 
„Ен 1:9 + е 9 = „в 2428. — гипергеометрическая функция, которую 
2 т 


согласно [6] можно для данных условий выразить через гамма-функции: 
= —1 
‚ЕР Ва + ‚а +2 - Е 2422] = Г(@ + 2 (1 гб ни(3)| . (25) 
т 
Если подставить (22), (24), (25) в (26) и произвести несложные преобразования, 
то можно получить окончательное выражение для интеграла [: 


—1 


ее А м _ @5 


Интегралы [,,[; и [. можно найти аналогичным образом. Окончательные их 
значения выражаются следующими соотношениями: 


рее Ца +1) мы" [> [я 


В 2В 


=] 


Ь= т = 1. (27) 


Подставляя зависимости (27) в (18) и производя несложные преобразования, 
получим выражение для искомого корреляционного момента: 


[+] 
5 0 +— 
г (28) 


й 
ь т 1 ит 
бои) --в[ Г >) ВГ (а) в 


Отрицательный знак корреляционного момента К (т, г/х) объясняется тем, что с 
повышением меры точности г рассеивание случайной величины т должно умень- 
шаться, а при малой мере точности г разброс параметров т должен быть большим. 

Если подставить значения параметров Т,,,,В и И,, определяемых выра- 
жениями (12), в (28), то можно определить апостериорное условное значение корре- 
ляционного момента К (т, г/ х) в зависимости от числа проведенных испытаний. 


Полученное выражение (28) определяет корреляционную связь случайных зна- 
чений математического ожидания и меры точности для нормального распределения 
статистических данных, которая должна учитываться при анализе и оценках стати- 
стических испытаний аппаратуры, что позволяет повысить их эффективность. 
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Если подставить значения величин Р(т / й. Р(+ 7 в К (т, г/ х) определяемых фор- 
мулами (14), (15) и (28) с учетом (12) соответственно в выражение (11), то можно 
найти зависимость функции условного апостериорного риска Р(ии,п2) от объемов 
испытаний и И› при заданных исходных значениях параметров ил, а и В ‚, опре- 
деляемых формулами (2), а также значений выборочных дисперсий 5 и $,, оп- 
ределяемых зависимостями: 

10 


1 = >), 52= "У н). 


пы И 
В этой сложной зависимости можно наблюдать следующие закономерности: 
1. При больших априорных значениях Б(т) и Г(ь) и малых (’, минимум риска 


соответствует большим объемам выборок и! и и.. 
2. При больших (’, и малых априорных значениях Б(м) и Б(ь) оптимальные 


величины и! и и› будут небольшими. 


3. Соотношение между оптимальными значениями 7, и И, существенно зависит 
от соотношений между выборочными дисперсиями 5 и 5, и коэффициентами потерь 
С;, Со. 

Для определения оптимальных значений и, и и> применение классического ме- 


тода приводит к сложным уравнениям больших степеней, решение которых возможно 
только численными методами. Это существенно усложняет анализ зависимостей опти- 
мальных значений И; от исходных данных. В данном случае определение оптимальных 


значений представляется целесообразным осуществить градиентным методом, сущность 
которого сводится к следующему: 
1. Выбирается начальная точка Мо. на поверхности функций р (и, у пз) ‚ соответст- 


вующая начальным значениям объемов выборок т ы п. й 
2. Определяются приращения аргументов: 
ь др(ин› из) 
д т 


ве др(пи из) 


Аи =- И2- 
д и> 


) 


где а — коэффициент пропорциональности. 
3. Определяются координаты новой точки Му: 
ей р+! _ РР 
т =Ш+Ат 12 =И2+Ап», 
где р -— номер шага итерации. 
4. Определяются условия окончания поиска: 


В др(пи.пь) др (ии, п) 
ди, ди> 
Например, если & > 01, то а =10 ‚ если 0,01 < # < 0,1, то а =1, если же = < 0,01, 


Е 


то поиск закончен. 
5. Проверяется достаточное условие существования минимума 


2 2 
д Р(ии*пь) > д р >пь) > 0. 


— ) 


Эт ди> 
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Расчеты показывают, что корреляционный момент К (т, г/ х) при и.>1 очень 


мал и его можно не учитывать, частные производные определяются при этом следую- 
щими зависимостями: 


др(п, п.) _ 5. Е - Е.А, 2Р,-25.В. 


ди Е! Е} 
др (п.п) _ 521 ЕзА 2Е.-25,В 
ди> Е Е 
где А=т+2а -2+2М№ 


В =4а +2М№ 
Ри =т(2а -2)+ Ме +2 —2)+ № 
Е. =2В +5 ии+ 5 п> 
№М=и:+п2. 
Вторые частные производные соответственно определяются соотношениями: 
д’Р(и, из) _ 2 Е! Е.- Аб. Е. + АЕ» _ 245: Е›-351 В | 


дп Е: Е} 
д’р (т.п) _ о Е:Е.-А5,Р.+А’Е» 2452Е._352В. 
д’п› Е! Е 


Таким образом, оптимальные значения объемов выборок и, и и› можно полу- 


чить с помощью приведенного выше алгоритма. 


Выводы 


1. При условии равенства выборочных дисперсий 5, и коэффициентов потерь 


С, целевая функция Р(и, п.) будет симметричной относительно координатных осей 
* 
Ой, и Оп» ‚ поэтому оптимальный объем выборки И, при первой серии испытаний 


будет равен оптимальному объему и › при второй серии испытаний. 
2. Изменения выборочной дисперсии 5, относительно 5, или коэффициента 
С; относительно (`. в 4 раза при прочих равных условиях вызывают изменение соот- 
ношения между оптимальными объемами выборок и! и и› не более чем на 15% — 20%. 
3. Увеличение коэффициента потерь (`, в 4 раза вызывает уменьшение значения 
оптимальной величины объема и, примерно на 10%. 


4. Уменьшение выборочной дисперсии 5, в 4 раза вызывает увеличение опти- 


* 
мального объема у; не более чем на 20%. 
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А 3 


А1-Аттоип Ай, А1-Аттоип Набап 
Тесйтсу ют Еуашаноп орше Оритит 51е оф\ай5иса 
Тез15 Бу Вауеяап Мейо45 


Опе оЁ Ше та]ог {ез{Ё гедитетеп юг ау оЁГ едшртепе Фа шсгеазез Фет 
еИслепсу 1$ ю тедисе Фе питбег оЁ {е5{5 аЁа слуеп ассигасу оЁ ай$Иса| 1пРогтавйоп. 
Вазе4 оп фе Вауезлап тефо4д, \Шсь {аКез шо ассоипЕ а рмом даа оЁ байзиса| за41ез, 
Фе а!оог т Гог езитаез оЁ Фе гези $ Гог банз$Иса[ (е5$ Ваз Бееп таде. 'ТЬ1$ а1еогифит 
ргоу14ез Ше гедге4 ассигасу \ и 1ез5 патЬег о (е55. 

Тре тефо4 Гог еуашаноп оЁ уашез Гог орита[ {езипе, \усьВ {аКез шо ассоцие фе 
шЙчепсе оф зеуега| Рас{ог$ фай гедисе Фе питбег о гедиге4 {е5$, Ваз Бееп оМашед Бу 
шшитпииайоп оЁ Вауезлап т1зК, свагацеттеа Бу Фе сВапое ш Фе уапапсе$ оЁ гапдот 
рагатеег$ оЁ ип4ег $ваЧу едитртеп ити» Фе {ез@п® ргосезз. Зиосе$е4 гези$ оЁ 
увез м\Ш ппргоуе Ше еНесНуепез$ оЁ %“айзйса1 {е3{5 оЁ Ше едитртепе. 

А$ а гези[ оЁ гезеагсй, ( \уаз Аеуе]оре4 ап а!еогит Гог еуашайоп оЁ Кеу рагатаегз 
Гог Ше (езйпе едитртепе ассог4те © Фе зайзисз$ Базе оп Ше Вауезлап тефоа, Ша 
$1осапйу гедисез Фе гедште4 питабег 0 {е5{5 аЁа слуеп ассигасу оЁ $анзиса| пРоптавоп. 

Тре тафетайса! ганоп$Ыр, \уысВ аейпез Фе сотеаНоп оЁ гапдот уааез оЁ те 
тафетайса! ехресайоп ап теазигез оЁ Ше ассигасу оЁ ромепом ФфзелБиНоп оЁ 
занзИса| ааа Гог Фе погта! 1а\у, Баз Бееп оМашед. 

Тре тешфо4 оЁ еуашайоп Гог оршта!| {е565 \аз 4еуе]оред. [ 15 Базе оп фе тит- 
пимабоп оЁ Вауезлап 115К ипаег диаагайс 1055$ Гапсйоп, \иеВ сВагасептез фе демлайоп оЁ 
Фе гап4от уааез юг ше шуезисае4 рагатаег гот 15 ауегасе уаше. ТБе заосе$е4 пле(- 
Вод аПо\у$ ‘акте ш® ассоипЕ фе шЙчепсе оф зеуега| ппроцапе Гасюгз оп фе оритла! $17е оЁ 
(2565 саглед оц т смо БасВез. 
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